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1 Introduction and notation 



A variety of important observables are essentially given by the correlators of currents 
with different tensorial structure. The cross section for e~^e~ annihilation into hadrons is 
governed by the vector current correlator, the Z decay rate by a combination of vector and 
axial correlators [Q]. Higgs decays can be expressed by the correlators of the corresponding 
scalar or pseudo-scalar current densities. Theoretical predictions for all these quantities 
in one- and two-loop approximation, corresponding to the Born level and the order a.^ 
corrections, are available since long. However, in view of the present or the foreseeable 
experimental precision improved calculations of these quantities are required, at least 
up to order a^, if possible even up to order a^. For massless quarks the NNLO result is 
available both for the vector p[ and the scalar correlators 0. Quark mass effects are often 
incorporated by inclusion of the lowest one or two terms in an expansion in m^/s. This 
approach is well justified for many applications f^. Even for relatively low energy values, 
down to about three to four times the quark mass this approach leads to an adequate result 
if sufficiently many terms are included in the expansion [^, |^ . Nevertheless it is desirable 
to calculate the correlator for arbitrary values of m^/s without directly invoking the high 
momentum expansion. Recently this was achieved for both real and imaginary parts of 
the three-loop polarization function with a heavy quark coupled to an external current. 
In a first step the method was applied to the case of the vector current correlator. In this 
paper also the axial- vector, scalar and pseudo-scalar cases are considered. This completes 
the relevant corrections to polarization functions for neutral gauge bosons induced 

by heavy quarks, more specifically to their "non-singlet" parts. Contributions from the 
double-triangle diagram, giving rise to "singlet" contributions, are not considered in this 
work and will be treated elsewhere. 

It is useful to define dimensionless variables: 

where q is the external momentum of the polarization function and s is the center of 
mass energy in the process e~^e~ — >■ hadrons. Then the velocity, v, of one of the produced 
quarks reads 

V = Vl-x^. (2) 

Every time the generic index S appears without further explanation it is understood that 
S represents one of the letters a,v,s or p. 

The polarization functions for the four cases of interest are defined by 

{-q'9,u + q,qu) U\q') + q.q^UHq') = t J rfxe^«^(0|T4(x)jf (0)|0) for 5 = v,a,{3) 

q^U\q^) = i J dxe''^^0\Tj\x)f{0)\0) for 6 = s,p, (4) 

with the currents 

fa = "^Tm^^ Ja = ^7m75''/', f = ^V', / = i^75^- (5) 
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In Eqs. and (^) two powers of q are factored out in order to end up with dimen- 
sionless quantities 11'' (g^). As we are only interested in the imaginary part the overall 
renormalization can be performed in such a way that this is possible. Furthermore we 
ensure that n''(0) = 0. The transformation from this scheme to other (overall) renor- 
malization conditions is discussed in App. ^ Concerning the renormalization it should 
be mentioned that for the scalar and pseudo-scalar current the combinations mj'^ and 
mj^, where m is the pole mass, have to be considered in order to arrive at finite results. 
Note that 11^ = and II^ is trivially obtained from 11^ through the axial Ward identity 
{q^fni{q^) = Am^q^ {W{q^) - q\dW{q^)/dq^)\g2=o)- 
The physical observable R{s) is related to n(g^) by 

R\s) = 12n ImU^q"^ = s + ie) ioi5 = v,a, (6) 

R\s) = Stt ImU^q"^ = s + ie) for5 = s,p. (7) 

It is convenient to define 

U\q') = nW'^(g2) + ^ii^C^n(i)'^(g2) + ( ^AlHY ui^l^^q^^ + . . . , (8) 

TT y TT J 

Yl(2),5 ^ dnj^'^ + CaCfU^^Y + CpTriiUP'^ + CpTUP'^ + CpTUf'^ (9) 

and similarly for R^{s). The abelian contribution 11^''''^ is already present in (quenched) 
QED and Hj^J originates from the non-abelian structure specific for QCD. The polariza- 
tion functions containing a second massless or massive quark loop are denoted by np"*'^ 
and n^^'^, respectively. U^'^ represents the double-triangle contribution. Our procedure 
will be applied to the first three terms in Eq. (^, the last two terms will be studied 
elsewhere. 

The paper is organized as follows: In the next section the expressions for the po- 
larization functions in the different kinematical regions are provided. In Section |^ the 
approximation method is described and the results are given in Section ^. The conclu- 
sions are finally presented in Section ^ 

2 Discussion of the kinematical regions 

This section provides a discussion of three kinematical regions where analytical results 
are available and contains the input data required for the approximation method. 

High energy region 

The input from the kinematical region where — ^ puts stringent constraints 
on the form of the polarization function and plays an important role for our procedure. 
In the limit of large external momentum the polarization function can be cast into the 
following form: 

n>0 ^ 
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The coefficients contain ln(— z)-terms up to third order. For the subsequental discus- 
sion only the terms with n = and 1 are needed. These ffist two terms can be calculated 
by simply Taylor expanding the diagrams in the mass Q. This leads to massless three-loop 
integrals for which the technique is known since long 0. The following decomposition of 
the coefficients is adopted: 



TT y TT J 

D^}f = E4'ilfln^) J e {0,1,2}, xe{A,NAl,F}, (11) 

where d^Jlff^ are numerical constants. For -0^°'*'^ and D!^^'^ the sum runs only up to A; = 1 
and k = 2, respectively. The results for the four cases are available in the literature 
0, H, ^. We have independently repeated this calculation. The results in the on-shell 
scheme are listed in Tab. Due to our renormalization condition (n'^(O) = 0) for a 
comparison with 0, §|, ^ the terms presented in App. ^ have to be taken into account. 
The analytic expressions for D"^ can be found in and for D'^ and in |^ . Dq and Df 
are listed in App. |B[ 

This information will serve as input for the procedure described in Section ^. It 
should be noted that in contrast to the vector case treated in p the axial-vector, scalar 
and pseudo-scalar correlator also develop cubic logarithms \n^{—z). If we had adopted 
the MS definition for the mass and fi^ = q"^, these cubic logarithms would vanish. 

Behaviour at = 

An important input to the behaviour of the polarization function originates from 
the Taylor expansion around = 0. In this case the three-loop diagrams have to be 
expanded in the external momentum leading to massive tadpole integrals. The calculation 
is performed with the help of the algebraic program MATAD written in FORM 0. It 
automatically expands in q up to the desired order, performs the traces and applies 
recurrence relations |Ty] to reduce the many different diagrams to a small set of master 
integrals. The structure of n(g^) is as follows: 



U\q') 






TT 


/^(2),5 
n 






+ ... , 



^Recently terms up to order n = 4 for the scalar and pseudo-scalar correlator and up to n = 6 
for the vector correlator have been calculated. This requires also the knowledge of massive tadpole 
integrals. 
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n = 


: 






n = 


-1 

: 1 




k 





1 


2 


3 





1 


2 


3 


^NA.n 

l,n 
0(2)." 


2.2222 
-3.9749 

1.3075 
-9.5651 

2.9723 
-1.2583 


-1.3333 
-1.0000 
0.1250 
-0.7175 
0.2306 
0.2306 


0.0000 
0.0000 
0.0000 
0.4583 
-0.1667 
-0.1667 


0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


2.0000 
0.0000 
2.7727 

-3.9903 
2.2899 

-5.1048 


0.0000 
-3.0000 
-0.3750 
-7.7083 
2.1667 
2.1667 


0.0000 
0.0000 
2.2500 
1.3750 
-0.5000 
-0.5000 


0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


„(0),a 

n(2)^i 

^A,n 

r,{2),a 

^NA.n 
n(2).a 


3.5556 
-2.0860 
-1.0387 
-9.8007 

3.6266 
-0.2639 


-1.3333 
-1.0000 
0.1250 
-0.7175 
0.2306 
0.2306 


0.0000 
0.0000 
0.0000 
0.4583 
-0.1667 
-0.1667 


0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


-2.0000 
1.9623 
-8.7928 
14.8142 
-8.5657 
6.5500 


2.0000 
1.5000 
0.2672 

-1.5224 
2.3606 

-4.5090 


0.0000 
-1.5000 
-1.3125 
-4.5417 
1.3333 
1.3333 


0.0000 
0.0000 
0.7500 
0.4583 
-0.1667 
-0.1667 


r){o).s 

^A.n 

n(2).^ 

^NA,n 
Din 

ni2):s 

^F,n 


5.3333 
-3.9623 
2.0416 
-28.2865 
15.5016 
-10.1876 


-2.0000 
-4.5000 

0.8578 
-5.2693 
-0.5273 

6.3423 


0.0000 
1.5000 
3.5625 
5.9167 
-1.8333 
-1.8333 


0.0000 
0.0000 
-0.7500 
-0.4583 
0.1667 
0.1667 


-3.0000 
7.8185 
-6.0543 
34.4008 
-21.3855 
18.5432 


3.0000 
3.0000 

-10.0450 
-7.3623 
9.3091 

-14.2997 


0.0000 
-4.5000 
-5.0625 
-12.9375 
3.7500 
3.7500 


0.0000 
0.0000 
4.5000 
1.3750 
-0.5000 
-0.5000 


n(o).p 

^ A.n 
^NA,n 

n(2)-p 
r,{2),P 

^F,n 


4.0000 
-2.9623 
1.0691 
-17.6658 
10.9195 
-11.5909 


-2.0000 
-4.5000 

0.8578 
-5.2693 
-0.5273 

6.3423 


0.0000 
1.5000 
3.5625 
5.9167 
-1.8333 
-1.8333 


0.0000 
0.0000 
-0.7500 
-0.4583 
0.1667 
0.1667 


1.0000 

2.6062 
-5.6088 

6.0218 
-0.9696 

0.8079 


1.0000 
-3.0000 
-6.8483 
-13.6485 

6.3253 
-3.5443 


0.0000 
-1.5000 

4.3125 
-2.4792 

0.5833 

0.5833 


0.0000 
0.0000 
1.5000 
0.4583 
-0.1667 
-0.1667 



Table 1: Numerical values for the coefficients Df^ at one-, two- and three- loop level for 
/i^ = m^. 

Although the calculation is performed analytically the results are listed in numerical form 
in Tab. |]with the choice /i^ = m^. The analytic expressions are given in App. |C[ 

For the vector correlator the first seven coefficients are already listed in 0, whereas 
all other results are new. 

Threshold behaviour 

At threshold it is most convenient to consider ffist the information about R^{v) and 
transform this subsequently into the corresponding expression for n'^(g^) via dispersion 
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(-,(2). 6 




i.ii 


^ 1 .11 


V 


1 

2 
3 
4 
5 
6 
7 

o 



1.06667 
0.45714 
0.27090 
0.18470 
0.13640 
0.10609 
0.08558 
0.0/094 


4.04938 
2.66074 
2.01494 
1.62997 
1.37194 
1.18616 
1.04568 

A no c c o 
0.93558 


5.07543 
6.39333 
6.68902 
6.68456 
6.57434 
6.42606 
6.26672 

0.10789 


7.09759 
6.31108 
5.39768 
4.69907 
4.16490 
3.74591 
3.40886 


-2.33896 
-2.17395 
-1.89566 
-1.67089 
-1.49436 
-1.35348 
-1.23871 
-1.14341 


0.72704 
0.26711 
0.14989 
0.09947 
0.07230 
0.05566 
0.04459 

O.OODYY 


a 


1 
2 
3 
4 
5 
6 
7 




0.53333 
0.15238 
0.06772 
0.03694 
0.02273 
0.01516 
0.01070 
U.UU /so 


1.70123 
0.71577 
0.39678 
0.25243 
0.17477 
0.12819 
0.09805 

A A 77 /I Q 


2.31402 

1.64824 
1.18988 
0.89978 
0.70769 
0.57409 
0.47720 

U.4U44y 


2.71368 

1.69284 
1.09011 
0.75793 
0.55831 
0.42925 
0.34095 

A 077CA 

U.z / /oU 


-1.04643 

-0.66383 
-0.43131 
-0.30145 
-0.22280 
-0.17167 
-0.13656 

A 1 1 1 

-U.iiioo 


0.34325 

0.06509 
0.02238 
0.01013 
0.00540 
0.00320 
0.00205 

A AA1 QA 


s 


1 

2 
3 
4 
5 
6 
7 
8 


0.80000 
0.22857 
0.10159 
0.05541 
0.03410 
0.02273 
0.01605 
0.01182 


0.45185 
0.77651 
0.52152 
0.35707 
0.25650 
0.19215 
0.14892 
0.11862 


0.03484 
1.42576 
1.52607 
1.32709 
1.11491 
0.93714 
0.79516 
0.68238 


-2.51105 
0.89546 
0.99889 
0.82750 
0.66459 
0.53778 
0.44174 
0.36853 


0.88148 
-0.35912 
-0.40023 
-0.33359 
-0.26913 
-0.21845 
-0.17982 
-0.15024 


0.71856 
0.19112 
0.07144 
0.03362 
0.01832 
0.01103 
0.00715 

C\ C\C\ A C){\ 

0.00489 


P 


1 
2 

3 
4 
5 
6 
7 
8 


1.33333 
0.53333 
0.30476 
0.20317 
0.14776 
0.11366 
0.09093 
0.07488 


2.33333 
2.61481 
2.12783 
1.75361 
1.48309 
1.28241 
1.12869 
1.00753 


2.71218 
7.03952 
8.27333 
8.46632 
8.33148 
8.09038 
7.82133 
7.55403 


-1.85805 
3.57843 
4.10409 
3.93931 
3.65904 
3.38071 
3.13004 
2.91007 


0.92593 
-1.23162 
-1.50248 
-1.47479 
-1.38655 
-1.29093 
-1.20154 
-1.12143 


1.31106 
0.49637 
0.25889 
0.16180 
0.11237 
0.08353 
0.06509 
0.05249 



Table 2: Numerical values for the coefficients at one-, two- and three-loop level for 

2 2 

— mr. 

relations. Whereas the treatment of the four cases in the high energy and small region 
is quite similar there is a big difference at threshold between the vector and pseudo-scalar 
correlators on one hand and the axial-vector and scalar correlators on the other hand. 
The latter are suppressed by a factor w.r.t. the former. This can already be seen by 
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considering the Born results: 



At 0{as) an expansion for — is helpful for the considerations below. It reads: 

= i?(0).-[^!^!li±^-4^+O(t;3), (14) 

iiW'-^ = i?(0).«|^!L!ll±^_2^ +0(t;5), (15) 

= /2W.^(^!llll±^_i^+0(t;5), (16) 



The exact results can be found in ITT] , 0. The analytical evaluation of the double-bubble 
diagrams with massless fermion loop insertions indicates that the characteristic scale 



of the first two terms proportional to vr^ is given by the relative momentum, the last term 
(which is due to hard transversal gluon exchange) by the mass of the heavy fermion. This 
motivates the decomposition adopted in Eqs. (p!^ - |T7D . 

Let us at Ola"^) first discuss the abelian contribution proportional to Cp. The ra- 
tio between R^J^'^ and R^^'^'^ is proportional to the Sommerfeld factor y/{l — e"^) with 
y = CpT^ds/v, which resums contributions of the form (tts/f)". Axial- vector and scalar 
correlators follow the P-wave scattering solution for the Coulomb potential. Hence an 
additional factor (1 + y"^ / {Air'^)) has to be taken into account in order to obtain the cor- 
rect leading term of 0{al) [1^. For the vector and pseudo-scalar contribution also the 
next-to-leading term in v can be determined by taking into account the correction factor 
arising from the exchange of transversal gluons which reads (1 — CpAas/T^) for the vector 
and (1 — C^Sas/Tr) for the pseudo-scalar case. As i?^"*'" and R^^^'^ already start at 0{v) 
the corresponding factors are not considered. Finally we arrive at 

RT- = 3 - 3.^) + flj'- = 3 (^^^^^.) + 0(v% (18) 

= 3 (^^^^«) + Oiv\ rT'' = 3 (i^ - |.^) + 0(v). (19) 

It is, of course, necessary to incorporate the strong 1/v singularity into the approximation 
method for the vector and pseudo-scalar case. In contrast the axial-vector and scalar 
current correlator are very smooth at threshold so that these terms are not used for our 
procedure. The comparison of these Fade results with the exact terms at threshold will be 
performed in Section ^ and demonstrates that the threshold behaviour is well reproduced 
— an independent test of our method. 
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The analytical results for the three-loop diagrams where a massless quark loop is 
inserted into the gluon propagator plus the corresponding real corrections are available 
for the vector current correlator and the remaining correlators as well[| in analytic 
form. The scalar case can also be found in 
leads to 



T6|. Expanding the results near threshold 



(2),a 



(2),. 



(2),P 



^(0),a 
i?(0),a 
^(0),P 



vr 



V 
V 
V 



In 



-In 

,6 





5 




18 




11 




18 




11 




18 








18 




(20) 



f211 



:22i 



(23) 



We include subleading terms proportional to In v in this expansion, since the agreement 
of our approximation improves visibly in those cases where the analytical result is known. 

In order to get the threshold behaviour for the non-abelian part it is either possible 
to use the QCD potential and the perturbative relation between Q;y(g^) and a;s(/i^) or to 



proceed as demonstrated in [T^ and deduce the gluonic double-bubble diagram, -R^^^'' 



from the corresponding fermionic contribution and evaluate it for a special choice of the 
gauge parameter ^. This is based on the observation that the terms proportional to Ca in 
the relation between a;y(g^) and ^^(/U^) are covered by the (one-loop) gluon propagator 
choosing ^ = 4. We will choose the second method since this trick is used also for the 
actual calculation. Following the expansion of the "double-bubble" result for ^ = 4 is 
taken to represent the expansion of the full non-abelian part. For the four correlators it 
is given by: 



R 
R 



NA 

(2),a 
NA 

(2),« 
NA 



-fin 



^NA 



^(0),s 









V 








V 








V 






V 



31 \ , 



11 



97' 



— In— + — +0{v 



24 



72 



11 v^s 97 , , 
-^ln- + ^)+0(.^ 



' 11 

r In — — 



24 



+ 



72^ 
31^ 
72 



0{v) 



(24) 
(25) 
(26) 
(27) 



To combine the results from different kinematical regions the above expressions for the 
imaginary part have to be transformed into analytical functions for Y\.^{q^) which respect 
Eqs. (H) and (|^). This can be done in close analogy to @]. 



^We would like to thank the authors of llSl for providing their results prior to publication. 
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3 The approximation procedure 



This section is devoted to the description of the approximation method. In order to 
save space we will not present explicit formulae. They look very similar to the ones for 
the vector case discussed in 0. The treatment of the abelian part of the pseudo-scalar 
correlator is in close analogy to [ |T8|] . 

In a first step a function Il'^(g^) is constructed which contains no high energy sin- 
gularities and no logarithmic terms at threshold. This is achieved with the help of the 
function 



G{z) = , u = ^ . 28 



The combination (1 — z)G{z) has a polynomial behaviour for z ^ and vanishes at 
threshold {z — > 1). For the case z — >■ —oo the expansion of G{z) develops logarithms 
starting with ln(— l/(4z))/(2z). This property is exploited and a function of the form 

J2 c„„^^^"(l-^)™(G(^))^ (29) 



n,ml 



where n is an integer and m, / > is constructed in order to remove the \n.{—z) terms 
of n^(g^). Whereas for the vector case described in no cubic logarithms appear (see 
Tab. 1^) and therefore quadratic combinations in G{z) are sufficient, for the other three 
cases this is not true: The axial- vector correlator develops \v?{—z)/z terms and com- 
binations like (1 — zY{G{z)Y are required. For the scalar and pseudo-scalar case also 
cubic logarithms appear which are not suppressed by powers of z, whence terms like 
z{l — z)'^{G{z)Y niust appear in the expression subtracted from W{q^). 

If logarithmic terms are present near threshold (nS'' and lif^'^) they are first sub- 
tracted and then the high energy singularities are removed. For the abelian polarization 
function where either 1/v singularities (n^"*'" and H^'*'^) or just constant terms (H^'*''^ 
and n^^-'''*) are present for z — > 1 only the high energy logarithms are removed. 

In a second step we perform a variable change. Via the conformal mapping 



= , ^=, z = = --. ( J(J) 

l + ^l-gV4m2 4m2 {1 + 

the complex g^-plane is mapped into the interior of the unit circle and the upper (lower) 
part of the cut starting at 2; = 1 is mapped onto the upper (lower) perimeter of the 
circle. The special points = 0,4m^,— cx) correspond to u = 0,1,-1, respectively. 
In this new variable we construct a function P{uj) for which the Fade approximation is 
performed. According to the different behaviour of 11'' (g^) near threshold actually two 
different functions have to be defined: 

P'i^) = (TT^^ - ' (31) 

P^V) = (ll(g^) - II(-oo)) . (32) 
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Thereby P^{uj) takes care of the cases where a 1/v singularity is present (H^^''' and 
■q(2),p'j^ The factor (1 — on) corresponds effectively to a multiplication with v. At this 
point we should mention that in order to incorporate also the constant term at threshold 
which has its origin in the correction factor introduced before Eq. (^) the combinations 
n^^^''' + 4n(i)''' and nj^'^ + 3n(i)'P are considered. The two-loop results H^^)'^ and H^^^'P 



may be found in |19[ and pOf, respectively. 

P^\ijj) treats all other cases where n^(g^) is just a constant for z = 1. Note, that 
this constant is unknown and consequently P^\l) may not be used for the construc- 
tion of the Pade approximation. P'^(l) is directly connected with the 1/v singularity 
and, of course, known. The high energy terms are treated in the same way for P^{uj) 
and P^\uj): Due to the subtraction of n''(— oo) the constant terms transform to -P(O) 
and the difference together with the prefactor 1/(1 + uoY projects out the 1/z sup- 
pressed terms in the limit uj ^ —1. Finally the moments from transform into 
derivatives of P{oj) at = 0. In total the following information is available for P\uj): 
{P^(-l), P^(0), P^'(i)(0), . . . , P^'(s)(0), P^(l)}. These eleven data points allow the con- 
struction of Pade approximations like [5/5], [6/4] or [4/6]. For P^^[uj) the threshold 
information P^^(l) is not available which means that at most Pade approximations like 
[5/4] or [4/5] may be constructed. 

For the non-abelian contributions proportional to CaCf there is an alternative ap- 
proach. Following the method outlined in the imaginary part of the gluonic 
double-bubble contributions, R^^^'^{s), can be computed analytically from the knowl- 
edge of the fermionic contribution Pp''''^(s). P^^)''^(s) is, of course, gauge dependent. 
However, for the special choice ^ = 4, where ^ is defined via the gluon propagator 
{— g'^" + ^ <f / (f) / {(f + ^e) the threshold behaviour of the non-abelian contribution 
(see Eqs. (p^ - pT]) ) and the leading high energy logarithms are covered by R^^^'^{s). 
Therefore it is promising to apply the procedure described above to the difference 
n^^'^(g^) — n^^^''^((3'^)|^=4 which has a less singular behaviour than n^^'^(g^). The re- 
sults for the non-abelian contribution presented in the next section are based on this 
method. 



4 Results 

After the Pade approximation is performed for the function P{uj) the corresponding equa- 
tions are inverted in order to get U^^q"^). In Figs. |l]-| the results are presented grouped 
according to the threshold behaviour. In Fig. |l| and ^ R{s) is plotted against the velocity 
for the vector and pseudo-scalar and the axial-vector and scalar case, respectively. The 
Figs. and ^ contain the corrections plotted versus x. Also the threshold and high energy 
approximations are shown (dashed lines). For the vector correlator terms up to (9(a;^^) 
are available Although in our procedure only terms up to (9(a;^) are incorporated the 
higher order terms are very well reproduced. For the scalar and pseudo-scalar polariza- 
tion function terms up to 0{x^) are available Again only the quadratic terms are 
build into the approximation method. However, the numerical coincidence with the high 
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energy approximations is very good. Actually it is hardly possible to detect a difference 
between the high energy terms and the Pade results when x is used as abscissa. A similar 



behaviour is observed for the axial- vector case where only quartic terms are available |21 



In this presentation it is not possible to notice any difference between the different 
Pade approximants. Minor differences can be seen after the leading terms at threshold 
are subtracted. This can be seen in Figs. ^ and It should be stressed that the vertical 
scale is expanded by up to a factor 100 in comparison with Figs. and ^. The following 
notation is adopted: All Pade approximations containing information up to Cg are plotted 
as a dashed line and the higher ones as full lines. The obvious exceptions are represented 
by a dash-dotted line and the exact results are drawn as dotted curves. 

The vector case is already discussed in 0. The inclusion of Cg into the analysis shows 
a further stabilization of the results. The plot for the abelian part in Fig. ^ contains 
altogether 14 Pade approximations. Eight of them contain information up to Cq (dashed 
lines), and six contain also information from C7 and Cg. These latter six lines coincide even 
on the expanded scale. The dash-dotted curve belongs to the [2/5] result and contains a 
pole for uj ^ 1.06. In the case of the non-abelian contribution 15 Pade approximations are 
plotted. Again the dashed lines belong to the lower order results. The dash-dotted curves 
are the results of two Pade approximants which have poles close to uj = 1 ([4/3] : uj ^ 1.07 
and [2/5] : uj ^ 1.06). 

For the pseudo-scalar correlator we find similar results concerning the behaviour of 
the Pade approximations when more information is included. The abelian contribution 
for the pseudo-scalar case contains 17 different results. The dash-dotted lines differ from 
the remaining ones significantly. The corresponding Pade approximants are [3/2] and 
[2/5]. A spread between the different Pade approximants can also be observed for the 
non-abelian contribution to the vector and pseudo-scalar cases. In both cases, however, 
convergence is visible if more information is included into the construction procedure. In 
the plot for 6R^^^J^{s), e.g., the dash-dotted line correspond to the Pade approximation 
[2/3] containing only the first three moments for — > 0. If this curve is ignored the spread 
is much less dramatic and the difference between the remaining Pade approximations 
shown is very tiny and completely negligible. The excellent agreement for the fermionic 
contribution in the pseudo-scalar case is comparable to the one for the vector correlator. 
In both cases the exact results , plotted as a dotted line, is indistinguishable from the 
approximations. 

Coming to the axial-vector and scalar case we would like to remind the reader that 
for these correlators no singularities are present in the limit v ^ 0. However, also here 
it is instructive to subtract the leading terms and look closer to the remainders SR^"^^'"" 
and 5R^^^''. As can be seen by comparing Fig. || and Fig. || the reduction in the scale 
lies between a factor two and ten. Also the very smooth behaviour of the subtracted 
results near threshold is clearly visible. One recognizes that, e.g. the remainders of the 
non-abelian and light-fermion contributions shown in Fig. ^are zero almost up to f ~ 0.2. 

Let us now consider the threshold behaviour of the abelian contribution. As mentioned 
in Section the corrections start with a term linear in v. We are now in the position 
to compare the Pade results with the exact expressions. In Tab. ^ the numerical values 
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P.A. 




P.A. 


A 






f3/3l 




19 /"^l 








[■^ /"^l 




f9 /4l 


O J- .4:4: 


[4/21 


29.69 


[3/41 


31.38 


[2/4] 


29.63 


[5/2] 


31.39 


[2/5] 


29.77 


[2/5] 


31.39 


[3/5] 


29.96 


[6/3] 


31.22 


[6/2] 


30.14 


[3/6] 


31.28 


[2/6] 


29.82 






[5/4] 


31.23 






[4/5] 


31.71 






[6/3] 


32.55 






exact 


31.75 


exact 


31.75 



Table 3: Comparison of the leading term at threshold for R]^ {6 = a, s) with the exact 
expression. All Fade approximants (P.A.) contain two terms from the high energy region 
as input. Only the number of moments from the expansion — > is different. 



of the coefficients of both the expansion and the exact result is shown. Although the 
analytically known terms are not incorporated into the approximation method they are 
very well reproduced by our method. 

In the abelian part of the scalar correlator there are two Pade approximants which 
differ significantly (dash-dotted lines) from the other eight results. One of them is a 
low-order Pade approximation containing only the information up to C2 and the other 
one ([4/3]) has a pole close to u; = 1 (1.02) which is reflected in the enhancement in the 
vicinity of the threshold. 

Both the non-abelian and fermionic contributions show an excellent agreement between 
the different Pade approximations. We should mention that at least 14 approximations 
are plotted and for i?P'''"(s) and rI'^^'^{s) in addition the exact results are included. Again 
no differences are visible. 

Finally we present handy approximation formulae for the abelian and non-abelian con- 
tributions. The procedure used to get them is described in 0. There the approximation 
formulae for the vector case are already listed. For completeness we repeat them at this 
point: 

,,,, _ (1^3^_,^a,. + ,H515_„32061 81 , 1-^ 

^ - o fi4 64 8 V y 1 + ^; 
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+ 100p=^/2(l - p) [2.21 Po(p) - 1.57 Pi (p) + 0.27 P2(p)] , (33) 
r.(2U,\ 351 o297 

^ l€=4 32 32 

+ 50p3/2(l -p) [1.73 Po(p) - 1.24 Pi (p) + 0.64 P2(p)] , (34) 

-^-§v + 18v' + ((27.; - 27.;^) (1 - ln2)) C(2) + (^ - v')) C(3) 

(1 - v^] In h iv - V^) 

32 V J J 1+v \ 16^ n l + v 

+ bQp^/^{l-p) [11.97 Po(p) - 25.37 Pi (p) + 20.67 P2(p) 

- 9.048 P3(p) + 1.85 P4(p)] , (35) 



+ 



16 ^ yy ' V 16 ^ V l + v 

+ 50p2/2^1 - p) [-1.88Po(p) + 3.31Pi(p) - 1.96P2(p) + 0.483P3(p)] , (36) 

1125 1779 o /189 261 , / ^ ^ ^\ 

-V + -7^^ + ( -^v - —v"^ + (-27?; + 45?;^) In 2 J ((2) 



64 64 V 4 4 

189 279 ,\ / 63 297 o\ , 1 - v 



) 3\ ^ / 63 297 2\ , 1 
-^)C(3) + (-- + -.^)ln- 



8 8 J ' V8 16 ; l + v 

( 243 297 o\ , 2 1 - 

+ V H v^ In 

V 16 16 y l + v 

2\ 3/2 



+ 50p=^/'(l -p) [1.30 Po(p) - 4.37 Pi (p) + 3.58 P2(p) - 0.91 P3(p)] , (37) 

.2^,1 99 147 o / 135 225 o /27 45 ,\ , A 

' +—V + v^ + v + v^ + [—V 'i;Mln2 C 2 

^ l?=4 32 32 V8 8 V2 2 ) ' 

/ 27 9 3\ / 45 135 2\ , 1 - 

+ -T?^ + TT:^' C(3) + -- r + In 



16 16 V4 8 ) l + v 

+ 50p3/2(l -p) [-3.94 Po(p) + 6.97 Pi (p) - 4.11 P2(p) + 1.00 P3(p)] , (38) 

il-v'^fT^^ ^ 2763 813 ,/ 9 63 , 

+^ 3PW'*' + V v^ +[ — V v^ 

V 4 64 64 V 4 4 

+ (9^ + 9^^) ln2) C(2) + (-ft; - fv') C(3) 

/81 63 2\ , 1 - ^ / 27 81 o\ , 2 
+ -r + T7:V In h -—v + —V^ In^ 



16 y l + v V 16 16 ; l + v 
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300 (^)- 



d(2),P 



100p^/2(l - j9) [2.79 Pob) - 1-83 Pi(j9) + 0.419 Psb)] , 



(39) 



459 

H - 

«=4 32 

27 9 o 

— V H 1;'^ 

16 16 



213 o 

1;"* - 

32 

C(3) + ft;^ln 



45 9 



2 

l-w 
1 + u 



2\ 3/2 



2v' 



50p3/^(l - p) [0.354 Po(p) - 0.251 Pi(p) + 0.456 P2(p)] , 



(40) 



where "p = 
polynoms: 



v)/{l + v), C(2) = vrV6, C(3) 



Po{p) 



1, Pl{p)=p, P2{p)- 

3 5q -^/\ 3 35^ 



;P + ^4(p) = ^ 



1.20206 and Pi{p) are the Legendre 

1 3 

2^ 

+ 



"2 + r • 



(41) 
(42) 



2" 2" ' 8 4 

^(2),'5 jg ^]-^g exact result for the gluonic double-bubble to be reconstructed from the 



fermionic contribution [IT? 



/?(2),5 



C=4 



-p(i)'^ 

3 



(43) 



For some cases the degree of the polynomial used for the fit has to be increased in or- 
der to end up with reasonable approximations. The first lines of the result contain the 
exactly known high energy and threshold contributions. The proceeding lines represent 
the numerically small remainder, R^^^'^'^^"^ with x E {A, NA}, which is plotted in Fig. ^ 
together with the result from the Fade approximation. 



5 Conclusions and summary 

The vacuum polarization function has been evaluated in order for vector, axial-vector, 
scalar and pseudo-scalar currents. The results take full account of the quark mass and 
are applicable between the production threshold and the high energy region. The method 
is based on the Fade approximation and uses the leading terms at high energies and at 
threshold plus the lowest eight coefficients of the Taylor series of n'^(g^) around zero. 
The stability of this approximation has been verified and excellent agreement between 
the present result and the predictions based on the high energy expansion is observed. 
These results can be used to evaluate the cross section for top pair production in electron 
positron annihilation through the vector and axial- vector current and the decay rate of a 
scalar or pseudo-scalar Higgs boson into top quarks in the full kinematical region and in 
next-to-leading order. 
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Appendix 

A MS definition of the polarization functions 

In this appendix we present the missing pieces needed to express the polarization func- 
tions, n'^(g^), in the MS scheme which means that in the expression obtained after the 
renormalization of as and m only the poles are subtracted. Expressing the results still in 
terms of the on-shell mass, m, 11'^ (g^) reads: 



(A.l) 



where the bar only refers to the overall renormalization. For the different cases we get 
{Cq is already listed in 0): 



Co 



vr V 4 



Cli^- -L^^ + (5 - 8 In 2)C(2) + ^((S) 



V144 8 



^ 14977 157 11 2 

+ ^^^^ 1^592" + 'm^''^ + 24^^™ 



;-- + 41n2)C(2) 



+ CpTni 



917 _ 14 _1 2 



^C(2) 



+ CpT [ -— - —L.^ - -L^ + -C(2) + -C(3) 



0, 



a, _ /67 
7r \36 
/131 _ 1 



~ + (5 - 8 In 2)C(2) + H^C(3)) 



+ - 

V TT 

+ CpCa + ]^L,m + ^^fim + {-\ + 41n2)C(2) 



127 
"96 



C(3) 



883 
576 



C(3) 



^ 72 27 



^C(2) 



C.T(-g-|L,„-iLj. + ^C(2)-^C(3) 



(A.2) 
(A.3) 



(A.4) 
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+ ( — ) Cp ^— - -i?4 + —L^rn - —Llr„ - -. 



, r3 

15 At™ 4 z*"^ 



+ (-y + 12 In 2) (1 - L^„) C(2) + (- 



+ CfCa (- 



1039 3 143 109 2 H 

, ~96~ + 4 " + 18" '^^ + ^^"^+12 

+ (2 - 6 In 2) (1 - C(2) + (f + ^L,^) C(3) - f C(4)) 
+ C^Tn, - ^L^^ - ^Lj^ - + (2 - 2L,^) ((2) + ^((3) 

+ CfT - ^L,^ - ^Ll^ - + (-4 + 4L,™) C(2) 

9 r / rnr o 



■12- ^L, 



27 



2^Mn.)C(3) + ^C(4); 



22 



3 



(A.5) 



+ 



505 3 13 15 2 3 3 

' KnJ f ^ V "6r^2 ^~T ^'"^16 '^"^^4 ^'" 

+ (I - 20 In 2) (1 - h,^) C(2) + (f + ^L,^) C(3) - f < 

/5429 3 33 109 2 II.3 

+ ^^^^ ['288 ' A ' Tq^"" ~ " 12 '^"^ 



■C(4)^ 
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g +2L^„)C(2) + -^ 



■C(4)j 

3^(^)] 



^,m) C(2) + -|C(3)] 



— 1 — 3L^m H — ^C'f i~n 
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+ 



. TT 



^ V 64 8 



+ (-^ + 61n2) (1 - 6L^„) C(2) - 27C(3) + yC(4)~ 

+ rr r 2641 3 163 251 2 ,11.3' 

27 \ 

+ (1 - 3 In 2) (1 - 6L^„) C(2) + 10C(3) - -((4) 

^ ^ ^161 11 

+ CpTni 

' V 48 2 



— T _ 1^7-2 _ 7-3 
2 ^iJim ^ -^yum 



.J„ + (1-6L^JC(2) + 4C(3)^ 
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/99 23 19 2 .3 
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4 Mm 
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TT 
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/ ^e; 7 4 1 A 
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15, 69 
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. .... 7 
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,2 /1409 
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15 _^r2 _3 3 



15 



^ 6 In 2) (1 + 2L^^) C(2) - 9C(3) + ^C(4) 



/ 129 1 ^ 185 251 2 11.3 
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+ (-1 + 3 ln2) (1 + 2L^„) C(2) + y C(3) - -C(4) 



/ 19 13 19 1 ^ 

+ Ci.Tn, ( - - - -L^^ - -Lj^ + (-1 - 2L^™) C(2) + -((3) 

1 



+ CfT[^^- fl,^ - ^Ll„, - + (2 + 4L,^) C(2) - yC(3) 



-|^2 

with Lf^rn = In/i^/m^. is Riemanns zeta-f unction with the values ({2) = 7r^/6, ({3) ~ 
1.20206, C(4) = 7r^/90, and -B4 ~ —1.76280 is a numerical constant typical for three-loop 
tadpole integrals pO[] . 
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B Analytic results for 

For completeness we present in this appendix the analytical results for Dq and Df. 
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-02 
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49 -0^ 
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3 4 /^^ y 

with ^(5) ~ 1.03693. appears in our results because of the normahzation condition 
n'^(O) = 0. 



C Analytic results for 

In this appendix we hst the first eight moments for 5^ — > expressed in terms of the 
on-shell mass, m, in analytic form for the four correlators. 
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For the vector case the first seven moments were aheady presented in 0. All other 
results are new. 
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Figure 1: R^"^^'"" and R^'^^'^ plotted against v. The dashed curves represent the threshold 
and the high energy approximations, respectively. Whereas for the vector case also the 
terms of order (w? / sY are available 0] for the pseudo-scalar correlator terms of order 
{w? / sY are plotted. 
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Figure 2: R^'^^''" and i?*^^)'** plotted against x. The dashed curves represent the high energy 
approximations including terms up to 0{x^'^) for the vector correlator and 0{x^) for the 
pseudo-scalar case. 
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Figure 3: R^'^^'"' and R^"^^'^ plotted against v. The dashed curves represent the threshold 
and the high energy approximations, respectively. For the axial-vector case terms of 
order (m^/s)^ are available and for the scalar correlator terms of order {m? / sY 1§] 
are plotted. 
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Figure 4: 

and plotted against x. The dashed curves represent the high energy 
approximations including terms up to 0{x^) for the axial- vector correlator and 0{x^) for 
the scalar case. 
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Figure 5: R^'^'^''" and i?^^)'? plotted against v. The leading threshold terms are subtracted. 
The dashed lines contain only information up to Cq whereas for the full curves also Cy 
and Cs is used. The obvious exceptions are represented by the dash-dotted curves. They 
are described in the text. 
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Figure 6: i?*^^)'" and R^'^^''^ plotted against v. The leading threshold terms are subtracted. 
The same notation as in Fig. |^ is adopted. 
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Figure 7: The remainder i?^^''"'''"^™, x G {A, N A) is plotted for different Fade approxi- 
mants (dashed hnes) together with the fit (sohd hne). 
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